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Abstract. We study the size of nodal sets of Laplacian eigenfunctions on compact Riemannian 
manifolds without boundary and recover the currently optimal lower bound by comparing the 
heat flow of the eigcnfunction with that of an artifically constructed diffusion process. The same 
method should apply to a number of other questions; for example, we prove a sharp result saying 
that a nodal domain cannot be entirely contained in a small neighbourhood of a 'reasonably 
flat' surface. We expect the arising concepts to have more connections to classical theory and 
pose some conjectures in that direction. 



1. Introduction 

We consider a compact n— dimensional C°°— manifold (M,g) without boundary. Writing A g for 
the Laplace-Beltrami operator, we are interested in Laplacian eigenfunctions 

— A g u = Ait. 

A natural object of study is the measure of the nodal set 

Z = {x € M : u(x) = 0} . 

An old conjecture of Yau [32] states that \Z\ ~ A'. For real-analytic (M,g) this was proven by 
Donnelly & Fefferman [7J, however, in the generic C°°— case even n = 2 is still open. In two 
dimensions, the best bounds are 

A^ < \Z\ < A*, 

where the lower bound is by Bruning [3] and Yau (unpublished), independently, and the upper 
bound is due to Donnelly & Fefferman [5] (with another proof given by Dong [5]). After exponen- 
tially decaying lower bounds were given by Hardt & Simon [TT] , the first polynomial bound is due 
to Sogge k Zelditch [IS] 

\Z\ > A 1 ^, 

although, as was later pointed out, an earlier result by Mangoubi |15j can be combined with the 
isoperimetric inequality to yield 

\Z\ > X^~^. 
Currently, the best bound is the following. 

Theorem. The volume of nodal sets satisfies 

U n - X ({x G M : u(x) = 0}) > A 2 ?. 

This was first proven by Colding & Minicozzi [5]. Subsequently, different proofs were given by 
Hezari & Sogge [11], Hezari & Wang [12] (for n < 5) and Sogge & Zelditch [19]. The arguments 
tend to be either local estimates on small balls in the style of Donnelly-Fefferman or global integral 
formulae. It is the purpose of this paper to give a new local approach exploiting the fact that a 
Laplacian eigenfunction behaves nicely under the heat flow. The approach is fully self-contained 
with the exception of our using a global inequality due to Sogge & Zelditch [T!5] 

l|Vw|| L oo(M) 

which is known to be sharp on spherical harmonics. We also sketch a variant of our proof that 
comes to rely on A= ^(m) si A~ ||u|| l°°{m) (also used by Sogge & Zelditch [H]), which is 
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easily seen to be equivalent because of ||Vu||ioo( M ) ~ A 1 / 2 ||u||itx 1 ( M ). 

As a by-product we show that, for c sufficiently small, a nodal domain cannot be contain in a 
cA -1 / 2 — neighbourhood of a sufficiently flat (n — 1)— dimensional surface in M. The statement 
is easily seen to be sharp because there are eigenfunctions on the flat torus T ra such that any 
nodal domains is contained in a CA -1 / 2 — neighbourhood of a geodesic of length f . As another 
by-product we formulate two very hard conjectures in the spirit of our approach whose resolution 
would imply a slightly sharper version of Yau's conjecture. 

2. Statement of results 

2.1. The main idea. The main idea is as follows: the heat equation with a Laplacian eigenfunc- 
tion as initial data and Dirichlet conditions on the nodal set has the explicit solution e~ u(x) 
(this is also the solution of the heat flow without any boundary conditions, however, we will be 
working locally). In particular, we have precise control on the rate of decay of the L 1 — norm in 
time. A natural candidate for comparison is the heat equation with the same initial data but 
Neumann conditions on the nodal set, which conserves L . However, the entire difference between 
Dirichlet and Neumann heat flow is caused by the nodal set and if it was too small, it couldn't 
account for the difference in behavior. 



Our proof is not actually using the Neumann solution because it requires some regularity on the 
boundary and would necessitate using reflected Brownian motion whose construction is nontrivial 
around the singular set. Instead, we choose another way and construct a stochastic process which 
might be interesting in itself: for small times it acts as a diffusion but as time grows the process 
converges back to initial data. We don't expect any serious obstacles if one were to work with 
actual Neumann solutions, some remarks on how to proceed are sketched in the end of the paper. 

Theorem 1. We have 

H' 1 - 1 ({x G M : u(x) = 0}) > A^. 

The (purely local) proof will give a sum consisting of local terms over all nodal domains D, the 
lower bound is then implied by a global inequality due to Sogge & Zelditch [19 

AV , NI , L1(D) > A„ IMI , L1(M) > X^. 
\\S u \\l-{d) ||Vu|| L oo (m) 

2.2. Geometry of nodal sets. We expect the argument to be applicable to more general diffu- 
sion processes and possibly other questions about Laplacian eigenfunctions. One example is the 
shape of nodal domains, where we briefly describe a simple geometrical result. It deals with the 
question whether nodal domains can be contained in a small neighbourhood of a 'flat' surface of 
codimension 1. In two dimensions, the statement reduces to a statement mentioned by Mangoubi 
[16) and to a theorem of Hayman [10] . however, in this generality it seems to be new. 

Let E C M be an arbitrary smooth (n— 1)— dimensional surface. We ask whether a nodal domain 
can be contained in a small neighbourhood of E. The e— neighbourhood of a generic geodesic 
(being itself as 'flat' as possible) on the torus T 2 already coincides with the entire torus - we 
therefore need to place some restrictions on E for the question to be meaningful. Using d g (-, •) to 
denote the geodesic distance, we call E admissible up to distance r if 

(1) Va; G M : d g (x, E) < r # {y G E : d(x, y) = d{x, E)} = 1 

(2) and, additionally, at any point x G E the surface can be locally written as the graph 
of a function <j> from the tangent plane in x to R with the spectral norm of the Hessian 
satisfying 

\\H4,\\ 2 <\. 

The first condition precludes the geodesies scenario while the second condition imposes local 
regularity on the surfaces. Both conditions could arguably be weakened, however, we are interested 
in the geometric implications of a 'smooth' theorem rather than optimal low regularity. 
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Theorem 2. There is a constant c > depending only on (M, g) such that i/EcM is admissible 
up to distance A -1 / 2 , then no nodal domain can be a subset of the cA -1 / 2 — neighbourhood o/S. 

As mentioned above, the function u(x) — Keexp(iy/\x) on T 2 endowed with the flat metric has 
all its nodal domains contained in a 0.5A -1 / 2 neighbourhood of a geodesic of length 1 (being ad- 
missible up to r = 0.5) and the example easily generalizes to higher dimensions. 

It is not difficult to see (from the proof) that the statement can be generalized to the union of 
admissible sets assuming they are sufficiently transversal at points of intersection and assuming 
the complement of the union does not have small connected components. Indeed, the proof 
immediately carries over to the following classical theorem (where the inradius is defined as the 
radius of the largest ball fully contained in the domain). 

Theorem (Hayman, [10]). There exists a constant c > 900 -1 such that for any simply connected 
domain Q C K 2 with inradius p 

P 

Note that the assumption of being simply connected is crucial but can be relaxed provided the 
'holes' aren't too small - this is the case for Laplacian eigenfunctions, where the lack of simple 
connectivity comes from other nodal domains, which cannot be too small themselves. 

2.3. Open questions. The main idea of getting information on the nodal set by constructing dif- 
fusion processes which deviate in behavior at the boundary - because its infinitesimally generated 
particles are absorbed/reflected differently - relies on having the right concepts. As it turns out, 
a rather crucial concept is the notion of heat content, which seems to have been studied by many 
people but in a predominantly asymptotic sense (see, e.g. [2]). We conjecture an isoperimetric 
principle for the heat content and show how such a principle would imply yet another proof of 
the main theorem. It also allows for a natural refinement of Yau's conjecture and seems to be the 
suitable for proving a Lieb-type generalization of Hayman's theorem - this is discussed in the last 
section of the paper. 

As for notation, the symbols < and ~ will always denote absolute constants depending only on 
the manifold (M,g). 

3. Proof of the main Theorem. 

3.1. Heat content. Given an open subset N C M, we use pt(x) to denote the solution to the 
following heat equation 

(d t - A g )p t {x) = xeN 
p t (x) = l xedN 
p (x)=0 xeN. 

The Feynman-Kac formula implies that this can be understood as the probability that a Brownian 
motion particle started in x will hit the boundary within t units of time. The quantity 

p t (x)dx 

is called the heat content of N at time t. It can be seen as a 'soft' measure of boundary size - 
for large times the function will be roughly of size 1 in the entire domain and all information on 
the size of the boundary will be lost. However, within t units of time a typical Brownian motion 
particle travels a distance of ~ t x l 2 . This can be immediately seen with Varadhan's large deviation 
formula [5T] 

Jim -4* log J£"(t, a;, y) = d(x,y) 2 , 

where K(-, •, •) is the heat kernel on the manifold (M, g). Hence, for very small times t, the function 
Pt(x) should therefore be mostly supported around the boundary and this yields a connection to 
the size of the boundary. 
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3.2. Definitions. Let D be an arbitrary nodal domain. Without loss of generality, we assume 
the eigenfunction u{x) to be positive within D: otherwise consider — u(x). Given u(x), we define 
a one-parameter functions v(t, x) as solutions to the heat equation with ^(a;)!^ as initial data and 
Dirichlet conditions. We set 

v{t,x) := e~ xt u{x) 

and note that v(t,x) then solves 

{d t - A g )v(t, x) = on D \ {u(x) = 0} 
v(t,x) =0 on {u(x) = 0} 
v(0, x) — u{x) on D. 

The Feynman-Kac formula for the Dirichlet problem is classical (see e.g. Taylor [ID])- Given an 
open domain ft € R n , f € C£°(ft), x € ft, then 

(e tAD f)(x) =E s (/(w(*))^(w.*)). 
where i > is arbitrary, denotes an element of the probability space of Brownian motions 
starting in x, E x is to be understood with regards to the measure of that probability space and 

Jl ifw([0,t])cfi 
ipn{ui, t) = < 

I otherwise. 

Here we see the connection with the heat content: for any point x € ft and any t > 

E a (Vto(w,t)) = l-p t (z). 

For reasons that will become apparent in the proof, for / £ Cfi°(fl) we define a second 'diffusion' 
operator 3 via 

(e tH /)(x) := E*(/(w(t))V>n(w, t)) + E x (l - ^(w, t))/(x). 
This operator is initially smoothing but ceases being so as time progresses. Indeed, if ft C M such 
that M \ ft is open, then if S is adapted to ft it is easy to see that for every x £ ft 

lim (e t3 f)(x) = f{x). 



e ta fdx = / fdx 



I p t (x)u(x)dx = / (1 - e tAo )u(x)dx. 
n Jn 



Finally, we claim that 
which is equivalent to 



A stochastic argument would be to say that among paths not leaving the domain, it is equally 
likely to start in a point x and end in a point y than the other way around - a statement that 
follows from the symmetry K(t, x, y) = K(t, y, x) of the heat kernel. 

3.3. A Comparison Lemma. We are interested in comparing the behavior of the Dirichlet 
solution e tAr> u with the behavior of e l "u on a fixed nodal domain D 1 where we assume without 
loss of generality that u\ „ > (otherwise: consider — u(x)). It is obvious from the definition that 

e ts u > e tA °u. 

Lemma. There exists a constant C > depending only on (M, g) such that 

e t3 u-e tAD u < Ct 1,2 Pt {x)\\Vu\\ L ^. 
Proof. The definition yields that 

e*"it — e tAD u = p t {x)u{x) . 

The quantity pt{x) is localized at a i 1 / 2 — neighbourhood of the origin and has an exponentially 
decaying tail at larger distance. The statement then follows from 

u(x) < d(sf,0D)||Vu|U«.. 

□ 
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3.4. Conclusion. Fix again an arbitrary nodal domain D and we assume again without loss of 
generality that ufsc)^ > 0. The heat equation, the comparison lemma and the L 1 — conservation 
of 2 give 

e~ xt / u(x)dx = / e tAn u(x)dx 

JD JD 

> / e ts u(x)dx - Ct 1/2 ||Vu|| L oo / p t (x)dx 

JD JD 



Therefore 



/ u{x)dx- Ci 1/2 ||Vu|| L ~ / Pt ( 

J D JD 

l-e~ xt \\u\\ L i (D) 



x)dx 



Pt(x)dx 



D 



t 1 ' 2 ||Vu|U» (2 ,) 



The result now follows from showing that up to constants depending on the manifold 



lim \ [ Pt (x)dx~H n - 1 (dD). 
t^o+ Vt Jd 



The involved quantities are essentially local and so is our argument: it is known that the critical 
set 

{x e D : u(x) = \Vu(x)\ = 0} 

has (n — 1)— dimensional Minkowski measure - recent results by Cheeger, Naber & Valtorta [3] 
even give bounds on its (n — 2)— dimensional Minkowski measure. Fix a small t > and cover the 
dD with cubes of side length t -1 / 2 . We call a cube regular if it does not contain an element of the 
singular set and singular otherwise. From the result above, it follows that the number of singular 
cubes is of order o(t^' n_1 - ) / 2 ^). 

Regular cubes. Let t > be fixed and let Q be a regular cube. Once a cube is regular, there is 
no need to further refine it as time tends to 0. To see this, we fix a second time-parameter z > 
and study the problem for z — > + locally in Q. Our desired result follows easily from Varadhan's 
large deviation formula [3T] 

lim — 4z \ogK(z, x, y) = d(x, y) 2 , 

2->0+ 

where K{-, •, •) is the heat kernel on the manifold (M,g). In particular, it implies that the bulk 
of the Green function is concentrated in a ~ z 1 / 2 — neighbourhood of the boundary dD n Q and 
using the regularity of the boundary, it is possible to compare hitting probabilities with transition 
probabilities from D to D c . This implies the statement. 

Singular cubes. It remains to show that the error introduced by those cubes containing an element 
of the singular set is small: it is not enough to note that their relative proportion is small because 
they are weighted with a factor t -1 / 2 , which becomes singular for small times. Using < pt(x) < 1 
gives 



^sing 



< ^l^singl 



Z} snl g consists of o(t^ 1 )/ 2 )) cubes of side-length i 1 / 2 , therefore 

^l^singl < ^ n/2 ^ ( "" 1)/2) ) =o(l). 

Improved estimates on the Minkowski dimension actually imply a faster rate of decay but these 
are not necessary for the conclusion of the argument. □ 
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4. Thin nodal sets 

In this section, we give a proof of Theorem 2. It is based on using the fact that at scale ~ A -1 / 2 the 
neighbourhood of a surface admissible up to A -1 / 2 behaves like the neighbourhood of a hyperplane 
in W l , which allows for problems to be reduced to well-known one-dimensional facts. 

Proof of Theorem 2. We consider the heat flow with Dirichlet conditions on the nodal domain D 

(d t - A g )v(t, x) = on D \ {u(x) = 0} 
v(t,x) =0 on {u(x) = 0} 
v(0, x) = u(x). 

The other case being identical, we assume u(x) to be positive in D. We start by proving 

Vi>0 inf p t (x) < 1 - e~ xt . 

We argue by contradiction to prove the following slightly stronger statement 

Vx E D u{x) = Pt(x) < 1 - e' xt . 

Given such a x € D, we see using the heat equation and Feynmann-Kac 

e~ Xt \\u\\ L ~ (D) = e- xt u(x) = E,(w(t)^(«(t))) 

< \H\L°-(D)^x{lpD(oj(t))) = \\u\\ L ao( D )(l -p t (x)). 

We now set the time to be t = A -1 . It remains to show that choosing c small enough, we can 
derive a contradiction to this bound on pt(x). Take a small c > and a cA~ 1//2 — neighbourhood 
of the admissible surface S. Assume D to be a nodal domain fully contained in that set and let 
x € D be such that u(ie) = HnH^o^m. The statement we need to contradict is pt(x) < 1 — e _1 
and we will do so but suitably bounding the probability of leaving the cA -1 / 2 — neighbourhood of 
S from below to achieve a contradiction. 

Since we are neither interested in the precise dependence of c on the constants in the regularity 
assumptions (which are set to 1 in the assumptions) nor in getting a good numerical value for c 
given these assumptions, it suffices to show that we are living at the right scale. From the regularity 
assumptions on E, we see that for a Brownian motion tu(t) starting in x the function dist(x,w(t)) 
behaves with respect to the probability of being contained in an interval on these scales comparable 
to a one-dimensional Brownian motion. Using B(t) to denote the 1-dimensional Brownian motion 
on R the relevant quantity is well- understood and the reflection principle (see e.g. [13) 1 implies 

P [ sup B{t) > cA~ 1/2 ) = 2P (b^X- 1 ) > c\- 1/2 ) . 

However, i?(A _1 ) is just a random variable following a normal distribution with mean fj, = 
variance a = A -1 / 2 , the contradiction follows from taking c sufficiently small. □ 

5. Comments and conjectures 

We believe the heat content to be a possibly valuable tool in the further study of Laplacian eigen- 
functions; this section studies some further implications, in particular we conjecture an isoperi- 
metric statement, which would imply another proof of our main theorem. 

5.1. Heat content isoperimetry. The heat content is a very stable notion and well-defined even 
for very rough domains. We consider the following statement to be highly plausible. 

Conjecture. Let (M,g) be a compact C°° — manifold without boundary. There exists a constant 
c > depending only on {M , g) such that for any open subset N C M and all times t > 

f p t (x)dx<<m. n - 1 (dN)Vi, 
Jn 

where the Hausdorff measure is understood to be oo if undefined. 
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Remarks. 

(1) Extremizers of the inequality need to have a smooth boundary: small irregularities in the 
boundary increase the surface measure but have very limited impact on the left-hand side. 
It would be interesting to understand the relation between the nature of extremizers and 
geometric properties of the manifold. Is there a connection to Cheeger sets? 

(2) If the domain N has the property that there is a real number r > such that each 
point x € N is contained in a ball of radius r (possibly centered around another point), 
then the two quantities should be comparable up to t ~ r 2 . If N is a nodal domain of 
the Laplacian, the Faber-Krahn inequality implies that the inradius is at most ~ A -1 / 2 
and therefore t ~ A -1 is the maximum time up to which we expect the quantities to be 
comparable. 

In particular, we conjecture that for a nodal domain both quantities are indeed comparable up to 
t = A . If this could be shown, it would immediately imply 

H n ~ 1 (x G M : u(x) = 0) < Ai 

5.2. Isoperimetry implies the main statement. Assuming the conjecture to be true, the 
second remark suggests that the maximum viable time for its application to a nodal domain 
without loss is given by t = A . We ignore possible issues arising in its construction (see, e.g. 
Bass & Hsu pQ) and assume the existence of the reflected Brownian motion on the nodal domain. 

Theorem 3. Assuming heat content isoperimetry and existence of reflected Brownian motion, we 
have 

({x e M : „(*) = 0}) > A* £ > A^\ 

where the sum ranges over all nodal domains D. 
We sum the local estimates (using [15] ) 

A* IM|l1(m) >A^". 

IMIl°°(m) ~ 

Proof. The proof has strong similarities to our previous argument. Again, without loss of gen- 
erality, we assume u(x) > on D and write e tA ° and e* Ajv for evolution under Dirichlet and 
Neumann data, respectively. Our new comparison estimate is even simpler and states that on a 
nodal domain D 

e tAN u-e tAD u< Pt (x)\\u\\ L - {D) . 

The proof for this comparison statement is easy to sketch: the difference between Dirichlet and 
Neumann solutions arises from those Brownian motions hitting the boundary. The difference is 
maximized if all those particles hitting the boundary arrive in a maximum of u after having been 
reflected. 



Integrating the comparison at time t — A 1 yields 

e _1 / u(x)dx — / e A °u(x)dx 
J d J D 



> I e N u(x) - p t (x)\\u\\ L ~>( D )dx 



u(x)dx — \\u\\l°°( D ) / p t (x)dx. 

D J D 

Therefore, assuming heat content isoperimetry and using the Sogge-Zelditch inequality 

n n -\dD)> x 1 / 2 -S^> A^\ 



□ 



STEFAN STEINERBERGER 



5.3. Yet another proof. Another variant of the proof is as follows. We assume again the exis- 
tence of reflected Brownian motion. The estimate 

e tA Nu _ e tA» u < Ct L ^pt(x)\\Vu\\ L o, 

follows from studying the difference between Brownian motion reflected and absorbed at the 
boundary and using the fact that within t units of time a Brownian motion may travel a distance 
of up to - i 1 / 2 . Then, 

e- xt [ u(x)dx = [ e xtAD u(x)dx 
Jd Jd 

> [ e xtAN u{x) - Ct l/2 p t (x)\\\7u\\ L . x{D) dx 
u{x)dx — Ct 1 / 2 !) V-uH^oc^) / pt{x)dx. 



D 



D 



In the limit t -> 0+ , 



hm — -= I pt(x)dx > At 



t^o+ ^/t J D ~ ||Vu||i-o 

and summing over all domains gives with the Sogge-Zelditch inequality the desired result. 

5.4. Geometric structure of nodal sets. The quantity pt (x) can be seen as a local measure of 
the closeness and size of the boundary. It seems extremely natural to conjecture the following. 

Conjecture. Let (M,g) be a compact C°° — manifold without boundary. There exists a constant 
c > depending only on (M,g) such that if p t (x) is globally defined with respect to the nodal set 
of a Laplacian eigenfunction with eigenvalue X, then 

P\-i(x) > c for all x £ M. 

This conjecture, while seeming likely, should be extremely difficult. In particular, combining it 
with heat content isoperimetry at time t = A -1 immediately gives 

H n - 1 {{x e M : u(x) = 0}) > Ai 

5.5. Heat content, Laplacian eigenvalues and the inradius. Given a domain £1 € W 1 , we 
can define the first eigenvalue of the domain as 

Al(Si )= i„ t &fi<4^ . 

feH£{n) J n f{x) 2 dx 

An inequality of the form 

where r is the inradius of the domain is trivial if n — 1 , true for simply connected domains in n = 2 
(Hayman's theorem) and false for n > 3. Indeed, in dimensions n > 3 it is possible to introduce 
very thin spikes making the inradius small but having little overall influence on the eigenvalue. 
However, Lieb [13] in a celebrated paper has shown that Hayman's theorem 'essentially' generalizes 
to higher dimensions: for any domain fi £ ffi™, there is a ball B of radius r ~ A -1 / 2 such that 
|f2 fl B\ ~ \B\ (the theorem also gives a precise relationship between the implicit constants). We 
conjecture that this phenomenon persists for the heat content. 

Conjecture. Let Q £ M™ be a bounded open set. If c\ > and some t > 

p t {x)dx > CiH^idtyVi, 



then there is a ball B of radius \ft such that \B fl f2| > <-a|.B|, where ci depends only on the 
dimension and c\ . 

Of course, via Lieb's theorem, this would establish a mutual equivalence between the first Laplacian 
eigenvalue, the size of balls having a large intersection with the domain and the time up to which 
heat content isoperimetry is sharp (t = Ai(il) _1 ). 
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